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Abstract 

In the Two-Higgs-Doublet Model(2HDM) with large extra dimensions(LED), we study the con- 
tributions of virtual Kaluza-Klein(KK) gravitons to 2HDM charged Higgs production, especially in 
the two important production processes e + e~ — ► H + H~ and e + e~ — ► H~~tb, at future linear collid- 
ers (LC). We find that KK graviton effects can significantly modify these total cross sections and 
also their differential cross sections compared to their respective 2HDM values and, therefore, can 
be used to probe the effective scale At up to several TeV. For example, at = 2 TeV, the cross 
sections for e + e~ — > H + H~ and e + e _ — > H~tb in the 2HDM are 7.4fb for m H - = 150 GeV and 
0.003fb for ran- = LI TeV and tan/3 = 40 , while in LED they are 12.1fb and O.Olfb, respectively, 
for At = 4Tev. 
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I. INTRODUCTION 



The idea that quantum gravity can appear at the TeV energyscale well below the Planck 
mass M p ~ 1.2 x 10 19 GeV was proposed in the 1990's . The large extra di- 

mensions(LED) modelQ] introduced by Arkani-Hamed, Dimopoulos and Davli has attracted 
much attention. It has been emphasized that the presence of large extra dimensions brings a 
new solution to the hierarchy problem, which can take the place of other mechanisms, for ex- 
ample, low-energy supersymmetry. However, it also interesting to examine a scenario which 
combines new physics beyond the Standard Model(SM) such as the 2HDm[t| and LED. 
This new possibility leads to different phenomenology than the usual LED scenario, which 
we explore here. 

In this extended LED scenario, as in the usual LED scheme, the total space-time has 
D = 4 + 5 dimensions. The SM and new particles beyond the SM live in the usual 3 + 
1— dimensional space, while gravity can propagate in the additional 5-dimensional space, 
which is assumed for simplicity to be compactified on the 5- dimensional torus T 5 with a 

e M p is no longer the relevant scale 



common radius R. Then the 4-dimensional Planck sea' 
but is related to the fundamental scale M„ as follows 



M; = M 6 s + \2tiR)\ (1) 
where M s ~ TeV. According to Eq. (fT)) . deviations from the usual Newtonian gravitational 



force law can be expected at distances smaller than R ~ 2.10 17 1CT* cm[8 



For 5 > 2, 



LED is consistent with the current experiments {9] since gravitational forces are not yet well 

probed at distances less than sbout a millimeter (However for 5 = 2, there are constraints 

fl 

arising from, e.g., supernova cooling, which require M s > 10 — lOOTeV if 5 = 2 [8]). 

2HDM LED can be tested at future high energy colliders. In 2HDM LED, just as in 
LED, there exist KK towers of massive spin-2 gravitons and scalars which can interact with 
the SM and beyond SM fields. There are two classes of effects that can probe LED: real 
graviton emission and virtual KK tower exchange. 

At future linear colliders the search for one or more Higgs bosons will be a central task. 
In the SM, the Higgs boson mass is a free parameter with an upper bound of mn < 600 - 
800 GeV (lOj. Beyond the SM, the 2HDM is of particular theoretical interest, in which the 
two complex Higgs doublets correspond to eight scalar states and spontanoues symmetry 
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breaking leads to five physical Higgs bosons: two neutral CP-even bosons h° and H°, one 
neutral CP-odd boson A , and two charged bosons H ± . The h° is the lightest and is a SM- 
like Higgs boson especially in the decoupling region (m^o 3> m z o). The other four are not 
SM-like and their discovery, particularly the charged Higgs bosons, would provide evidence 
for the 2HDM. If the H bosons have mass m H ± < m t — m&, they will be produced mainly 
through the t — > bH + decays of top quarks, which can be produced singly or in pairs at an 
e + e _ LC [llj|. If there is sufficient center of mass energy available, \/s > 2rriH±, then charged 
Higgs pair production, e + e~ — > H + H~, will be the dominant production mechanism 11]. 
However, if mw± > max(mt — mj, \/i/2), then H ± bosons can only be produced singly. And 
e + e _ — > //"~£ft|l2Lll3{| is one of the most important single charged Higgs production processes 
that can also be used to measure the relevant Yukawa couplings. Therefore, for this process, 
we will focus on the regime where y/s/2 < < y/s — m t — m&. For y/s = 500 GeV 

(1000 GeV), this implies that 250 < m H ± < 320 GeV (500 < m H ± < 820 GeV). 

There are several kinds of such 2HDM models. In the model called type I, one Higgs 
doublet provides masses for both the up-type and down-type quarks. In the type II model, 
one Higgs doublet gives masses to the up-type quarks and the other one to the down-type 
quarks. In the type III model, both doublets contribute to generate the masses for up-type 
and down-type quarks. In this paper, we will concentrate on the type II 2HDM(2HDM II), 
which are is favored by the Minimal Supersymmetric Standard Model (MSSM) 14 1. We will 
also briefly discussed the results in the type I 2HDM. 

In the following we consider the the contributions of virtual KK gravitons to the 2HDM 
charged Higgs production, especially in the two important production processes e + e~ — > 
H + H~ and e + e~ — > H~tb at future linear colliders. The presentation is organized as follows: 
In Sect. II we present the calculations. In Sect. Ill we give the numerical results and discuss 
them. Sec. IV contains a brief conclusion. 



II. ANALYTIC CALULATIONS 

In this section we derive the cross section for 2HDM charged Higgs production. The 
2HDM diagrams, including the additional virtual KK gravitons (G$), which contribute to 
the processes e + e~ — ► H + H~ and e + e~ — ► H~tb are presented in Figs. 1 and 2, respectively. 
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FIG. 1: The 2HDM Feynman diagrams and graviton mediated diagrams for e + e — > H + H at 
the tree level. 
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FIG. 2: The 2HDM Feynman diagrams and graviton mediated diagrams for e + e ^ H tb at the 
tree level. 

In the four dimensional description the interaction Lagrangian between the scattering 
fields and the KK gravitons (G^J) or KK scalars (H^) is given by [l^j ] 



int 



M, 



1 Wcgr^--J 3(n V 2! 



n + 2 



At > 



(2) 



where = (m, n 2 , .., r^) with n^'s being integers, M p = M p /\/&k ~ 2.4 x 10 18 GeV is the 
reduced four dimensional Planck scale, and is the energy-momentum tensor of the scat- 
tering fields. The ra-th KK mode graviton and scaler masses squared are both characterized 
by m?^) = \n\ 2 /R 2 . Since the trace of the energy-momentum tensor is proportional to the 
mass of the fields due to the field equations, we neglect processes mediated by the KK scalars 
in our study for the linear e + e~ collider, and consider only the processes mediated by the 
KK gravitons. 

From Eq .© we can derive the relevant Feynman rules to be used in our calculations, 



which can be found in Ref. 



unitary gaimeQ? ii> given by 
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161 ] . The numerator of the graviton propagator pv- va $ j n the 



/i/3 not 



+ ..., 



(3) 



where if" is the Minkowski metric. The dots represent terms proportional to the graviton 
momentum q^, and since q^T^ = 0, give a vanishing contribution to the amplitude. Note 
that the the numerator of the graviton propagator in Ref. 3| is twice Eq. (jHJ) and is the 
same as presented in Ref. However, the coupling constant squared in Ref. [l^J, i.e., 

(\/87r/Mp) 2 in Eq. © is twice the one in Ref. [16j ]. Hence the results in the two references 
are consistent and we have been careful of this issue in our calculations. 

We can write the amplitudes for the Feynman diagrams shown in Figs.l and 2 as follows: 

^2 ;„2(„2 n 2 



M 2a 



M 



26 



M 



2c 



^v( Pl )(fa 2 - V\)u{ P2 ) + f (cl f } 2 S(pi)(fl,- - sl)P L - s 2 w P R ]u(p 2 ) 

s 2{s - m 2 z )c 2 w s 2 w 2 
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—ie 2 g 



3\/2m w s[(ki + k 3 



mi 



m b tan f3P R + m t cot (3Pl 



L + fa + m b )^v(k 2 , m b ) + 



-ie 2 g 



V2m lu (s - mjjc^s^lih + k 3 ) 2 - m%\ 2 



m b tan j3P R + m t cot /3Pl 

-2ie 2 # 
3 v / 2m w s[(A; 2 + fc 3 ) 2 - m?] 

m b tan /3Pr + m t cot (5Pl 

ie 2 g 

y/2m w (s - mf j c 2 iu s 2 v [(k 2 + k 3 ) 2 - m 2 

1 2. n 2 
2- 5 4)ft+3 



h+ fa + m b ) 7 "[(^ - -si )P L - -^s 2 PflJvCfca, m b ) 



u(Pi)7/x«(p2)«(fci, mh' 1 (fa+ fa - m t ) 



v(k 2 ,m b ) + 



2 APi)l^ ~ s l) p L - 4 p fl]«(i ? 2)-«(A;i ) m t ) 
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^ - -sl)P L + ^s 2 w P R ]{fa+ fa - m t) 



-ie 2 g 
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-ie 2 g(c 2 w - s 2 ' 
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2c 
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M 



2/ 



4v / 2m U) [(fci + k 



ig gp^«P x d 

v{pi)C flu u{p 2 



m 



h- 



M 



2fl 



u{k 1: m t ) 




m b tan (3Pr + m t cot /^P^ 



v(k 2 ,m b ) 



(4) 



with 



-1 



E 



i 



Aa/3 



m p T s ~ m f. 

- + (/X <-> ^)] 

m 2 H -r] af3 - {^(k'zY {r]^ a r] uf 3 + r} VQ r}^ - r]^ v r} a p] 
[la{h + h- k 2 ) p - r] a p(^ 1 + k/ 3 - k/ 2 - 2m h )\ + (J3 «-> a) 
[7« (^i -h-k 2 )f3- VapWi- 2m t )\ + {f3 ^ a) 

m 2 H -r] af3 - £^(£4 + fc 2 ) w [77^^^ + T\ va f]^ - ri^Vap} , 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



where s = (pi+p 2 ) 2 , Pl = (1-75)/ 2 , Pr = (l+7 5 )/2, s w = sin w , c w = cos w and tan (3 is the 
ratio of the two vacuum expectation values in the 2HDM. Note that the contribution from 
Feynman diagram Fig. 2(g) vanishes since the trace of the graviton appears in this diagram 
and, therefore, in the limit of vanishing electron mass, this contribution also vanishes 17]. 

Q in Eq. (jSJ) represents the summation of the KK excitation propagators. If the summa- 
tion over the infinite tower of the KK modes is performed, one will encounter ultraviolet 
divergences. This happens because LED is an effective theory, which is only valid below 
an effective energy scale. In the following calculations we naively introduce an ultraviolet 



cutoff for the highest KK modes and replace the summation by 

1 



15] 



4tt -1 ^ 



(11) 



where is a cutoff scale naturally being of the order of the fundamental scale M s 
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Finally, the cross sections for the charged Higgs production processes in LED following 
from the amplitudes are: 



a(e + e 



> H + H-) = j-J d^\M la + M lb \ 2 (12) 
o(e + e- -> H~tb) = — j d^\M 2a + M 2b + M 2c + M 2d + M 2e + M 2/ + M 2 /, (13) 



and the final state phase space hypercube elements are defined as 



2 d 3 k , \ / 2 



3=1 



d »»=(n^) us) 

III. NUMERICAL RESULTS 

In the numerical calculations, we used the following set of SM parameters [3]: 

a ew {m w ) = 1/128, m w = 80.419GeV, m t = 178GeV, m z = 91.1882GeV. (16) 

For the Yukawa coupling of the bottom quark at the H~tb vertex, as suggested by Ref.jl^. 
we used the top quark pole mass and the QCD improved running mass rrib{Q) with m b (m b ) = 
4.25 GeV, which were evaluated using the NLO formula ^| as follows: 

m b (Q) = U e (Q, m t )U 5 (m t , m b )m b (m b ) , (17) 

Here the evolution factor Uf is 

\a s {Qi)s L 47r J 

d (f) = _JL_ jU) = 8982 - 504/ + 40/^ 

33-2/' 3(33 -2/) 2 ' 1 ' 

and / is the number of active light quarks. For the energy scale Q, we chose Q = 
as in Ref. [13j for e + e — > H tb. 

A. e+e- -» H+H- 

In Fig. 3 we show the dependence of the the cross sections in 2HDM II and LED for 
the process e + e~ — ► H + H~ on m H - assuming ^/s = 1000 GeV. The dashed lines represent 



the 2HDM II results and the solid lines represent the cross section in 2HDM II LED for 
different values of A^. From the figures one can see that when ran- is small(< 400 GeV) 
the LED effects can be very large up to A^ = 2 TeV. For example, when \fs = 1 TeV and 
m H - = 200 GeV the cross section in the 2HDM is 23.8fb, while for 2HDM II LED it is 
I52.4fb for A r = 1.5 TeV and 36.6fb for A T = 2 TeV. 

In Fig. 4, the cross section in the 2HDM II and LED for the process e + e~ — > H + H~ 
are plotted as functions of yfs assuming ra#- = 200 GeV. From this figure we see that 
the effects of virtual graviton exchange can substantially modify the e + e~ — » H + H~ cross 
section compared to its 2HDM II value, especially at large \/s(> 4 TeV), enabling LED up 
to At = 8 TeV to be probed. 

In Fig. 5 we show the dependence of the cross sections in the 2HDM II and LED for 
the process e + e~ — ► H + H~ on the effective energy scale Ay assuming ran- — 150 GeV and 
y/s = 1000 GeV, 2000 GeV, and 3000 GeV. As expected, when \/s is fixed the cross section 
in LED tends to the value in the 2HDM II when At tends becomes very large. In addition, 
the figure also shows that for large \/s and small At the LED effects are very large. For 
example, for TeV and At=4 TeV the cross section in the 2HDM II is 3.4fb while the 

value in LED is 58.7fb. 

In Fig. 6, we present the differential cross section da/dcosO in the 2HDM II and LED 
for the process e + e~ — > H + H as a function of cos#, where the scattering angle 9 is the 
angle between H~ and the incoming positron assuming ran- = 200 GeV, \fs = 1000 GeV, 
and At = 2000 GeV. Both differential cross sections in the 2HDM II and LED vanish 
at cosd = ±1. However, their shapes are different. Note that the LED cross section is 
not symmetric under cos 9 <-» — cos 9 while the 2HDM II one is symmetric. We further 
have plotted the different contributions to the LED cross section and find the asymmetry 
comes from the interference of the LED and 2HDM II amplitudes, which is in proportion 
to (t — u){tu — m H -)/A T { t and u are Mandelstam variables) and thus to cos#sin 2 # . 
Moreover, the LED contribution has peaks in both the forward and backward regions and 
also vanishes at three points: cos# = and ±1, as found in Ref. |20 ! ], since it is proportional 
to (t — u) 2 (tu — ra H -)/A T and thus to cos 2 9 sin 2 9. 
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B. e+e" -> iJ-tb 

In Fig. 7, we show the dependence of the the cross section in the 2HDM II and LED for 
the process e + e _ — ► H~tb on mg- assuming y'i = 1000 GeV for tan/3 = 40. The dashed 
lines represent the 2HDM II results and the solid lines represent the cross sections for LED 
with different values of A^. This figure shows that when is small(< 2.5 TeV) and m#- is 
also small(< 600 GeV), the LED effects are very large. For example, when tan/3 = 40 and 
m H _ = 550 GeV the cross section in the 2HDM II is 0.017fb. 

In Fig. 8, we present the dependence of the cross section in the 2HDM II and for LED 
for the process e + e _ — > H~tb on tan (3 assuming y/s = 1000 GeV and mjj- = 520 GeV. 
Observe that the cross section in both the 2HDM II and LED exhibit minima close to 



tan (3 ~ Jmt/mb ~ 6, since the average strength of the tbH coupling, which is proportional 



to yuif cot 2 (3 + m 2 tan 2 /3, is then minimal 2l|. From these figures one also sees that when 



A^ = 2 is small(< 2.5 TeV) and m H - is also small(< 600 GeV) the LED effects are large for 
= 1 TeV. 

In Fig. 9 the cross sections in the 2HDM II and LED for the process e + e~ — > H~tb 
are plotted as functions of y/s, assuming m H - = 800 GeV and tan/3 = 40. Since we are 
interested in the case where y/s/ 2 < m#± < y/s — m t — m^, y/s should be in the range from 
about 983 GeV to 1600 GeV. From this figure we see that, as in the case of e + e~ — ► H + H~ , 
the effects of virtual graviton exchange can significantly modify the e + e~ — > H~tb cross 
section compared to its 2HDM II value, especially at large \/s(> 1.5 TeV), where LED up 
to about At = 3.5 TeV can be probed. 

In Fig. 10 we show the dependence of the cross section in the 2HDM II and LED 
for the process e + e~ — > H~tb on the effective energy scale A^ assuming m#- = 
260, 550, andllOO GeV, and y/s = 500, 1000, and 2000 GeV, respectively. As expected, when 
y/s is fixed the LED cross section tends to the value in the 2HDM II when A^ becomes very 
large. This figure also shows that for large y/s and small At the LED effects can be very 
large. For example, when y / i=2TeV and Ar=4TeV the cross section in the 2HDM II is 
0.006fb, while in LED it is 0.013fb. 

In Fig. 11 we show the differential cross sections da / 'dPr(H~) in the 2HDM II and LED 
for the process e + e~ — > H~tb as functions of the transverse momentum of the charged Higgs 
boson P T assuming m H - = 520 GeV, ^ = 1000 GeV, and A T = 1500 GeV, 2000 GeV, 
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2500 GeV and 3000 GeV. The shape of the differential cross sections in the 2HDM II and for 
LED are slightly different. The LED effects generally enhance the differential cross sections, 
especially when is small. 

We briefly note the relevant results in type I 2HDM. The only change is in the fermionic 
couplings. Thus, the 2HDM and LED results for e + e~ — > H + H~ remain the same. However, 
the 2HDM results for e + e~ — > H~tb in the type I model are different compared with the 
type II model. For example, the results in type I model will decrease with increasing tan /3, 
since the H~tb coupling is proportional to cot (3 in the type I model. However, the LED 
effects will be similar: large effects when \fs is large and At is small. 

Finally, we briefly discuss the experimental sensitivity. We will assume that the LED 
effects can be observed, if j3] 

Ac- = a LED - <7 S M> £ , (19) 

where the integrated luminosity £ are assumed to be 500/6 -1 . For e + e~ — > H + H~, from 
Fig. 5, we see that the LED can be probed, for example, up to 2800 GeV and 7200 GeV with 
mu- = 150 GeV for y/s = 1000 GeV and 3000 GeV, respectively. For e + e _ — > H~tb, from 
Fig. 10, we see that the LED can be probed, for example, up to 1400 GeV and 3600 GeV for 
m H - = 260 and 1100 GeV, and ^ = 500, and 2000 GeV, respectively. 



IV. SUMMARY 



In the 2HDM II with large extra dimensions we investigated the contributions of virtual 
KK gravitons for the 2HDM II charged Higgs production in the two important production 
processes: e + e~ — > H + H~ and e + e~ — > H~tb at future linear colliders. We found that 
KK gravitons can significantly modify these total cross sections and their differential cross 
sections compared to the corresponding 2HDM II values and, therefore, can be used to probe 
the effective scale At up to several TeV. For example, at y/s = 2 TeV the cross sections for 
e+e" -> H+H and e + e~ -> H'tb in the 2HDM II are 7.4fb for m H - = 150 GeV and 0.006fb 
for mu- — 1.1 TeV with tan/3 = 40, while in LED they are 12. lib and 0.013fb, respectively, 
for A T = 4 TeV. 
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FIG. 3: Dependence of the cross section in the 2HDM II and LED for the process e + e — ► H + H 
on m H - , assuming y/s = 1000 GeV. 
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FIG. 4: Dependence of the cross section in the 2HDM II and LED for the process e + e — ► H + H 
on y/s, assuming m H - = 200 GeV. 




FIG. 5: Dependence of the cross section in the 2HDM II and LED for the process e + e — > H + H 
on the effective energy scale At, assuming m H - = 150 GeV, \fs = 1000 GeV, 2000 GeV, and 
3000 GeV. 
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FIG. 6: The differential cross section da/dcos6 in the 2HDM II and LED for the process e + e~ — > 
H + H- as functions of cos 9, assuming m H - = 200 GeV, y/s = 1000 GeV, A T = 2000 GeV. 




FIG. 7: Dependence of the cross section in the 2HDM II and LED for the process e + e — ► H tb 
on m H -, assuming y/s = 1000 GeV, tan/3 = 40. 
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FIG. 9: Dependence of the cross section in the 2HDM II and LED for the process e + e H tb 
on y/s, assuming m#- = 800 GeV, tan/3 = 40. 
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FIG. 10: Dependence of the cross section in the 2HDM II and LED for the process e + e — > H tb on 
the effective energy scale A T , assuming m H - = 260, 550, 1100 GeV, and y/s = 500, 1000, 2000 GeV, 
respectively. 
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FIG. 11: The differential cross section da/dPx(H^) in the 2HDM II and LED for the process 
e + e~ — > H~tb as a function of the transverse momentum of the charged Higgs boson Pj>, assuming 
m H - = 520 GeV, ^/s = 1000 GeV, A T = 1500 GeV, 2000 GeV, 2500 GeV and 3000 GeV. 
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